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ABSTRACT 
The use of Laplace Transforms has long been one of the primary 
methods utilized in solving linear differential equations but recently the 
extension and application of this method to the solution of basic non- 
linear systems has been proposed. This latter technique is presented 
and expanded to include various non-linear functions. The applicability 
of this method to systems of varied form and complexity is then explored 


and the validity of the derived solution investigated. 


DUDLEY KNOX LIBRARY 


NAVAL POSTGRADUATE 
SCHOO 
MONTEREY, CA 93943-5104 ; 


TABLE OF CONTENTS 


I. INTRODUCTION -w------terr ccc cc ccc ccc ctr ccrco 5 
Il, LAPLACE TRANSFORM CONCEPT -----2---2r-ceccrcrc rrr F 
ete MEO D Or INTEGRATING BY PARTS \=——=--=----- 6 
B, APPLICATION TO DERIVATIVES ------------------- 7 
em PPLiG2TION TO NON-LINEAR FUNCTIONS ----=----- cz, 
Ii, LAPLACE TRANSFORMS FOR NON-LINEAR FUNCTIONS ------ 9 


A, TRANSFORMS FOR THE PRODUCT OF A FUNCTION 
Aah oobi enna SSS See 9 


B, TRANSFORMS FOR THE PRODUCT OF A FUNCTION 
Aeon een DER AMVE ==> Ss —- = SSS yz 


Ge TRANSFORM OF THE PRODUCT OF POWERS OF A 
Ree hie AND Fo WEkiy Ai hae = > SS > SS SS SS = 14 


D. THE TRANSFORM OF THE FIRST DERIVATIVE RAISED 
ROMER OPOWER =a2eeo- 5225 S600 222 Seesc4 eS tke 


IV. APPLICATION OF LAPLACE TRANSFORMS TO NON-LINEAR 
SYSTEMS ----------------------------------- 17 


A, EXAMPLE OF A FIRST-ORDER NON-LINEAR SYSTEM, 
EQUATION 1 ------------------------------ 17 


B, EXAMPLE OF A NON-LINEAR EQUATION WITH SQUARE 


OF DERIVATIVE, EQUATION 2 ------------------- 20 

C. EXAMPLE OF THE PRODUCT OF THE VARIABLE AND ITS 
DERIVATIVE, EQUATION 3 ---------------------- 25 
V¥V, CONCLUSION -------------------------------- 29 
APPENDIX A COMPUTATION OF INVERSE LAPLACE TRANSFORM --- 31 
APPENDIX B LIST OF LAPLACE TRANSFORMS --------------- 35 
LIST OF REFERENCES ------------------------------ 36 


INITIAL DISTRIBUTION LIST === eee 


FORM DD 1473 





Sa 





¥ 
a 7 7 4 7 i 7 
—_— — os 
re ac a ze “al 


I, INTRODUCTION 

Laplace Transforms have been typically applicable and limited to 
linear systems. However, Baycura [Ref. 1 and 2]in two recent papers 
has proposed an extension of the Laplace Transform method to non-linear 
systems. By expanding the Laplace integral into an infinite series and 
integrating by parts n - times an infinite series was developed for non- 
linear terms. Through a further correlation of these derived series with 
series of linear transforms the transform for a non-linear term was obtain- 
ed. These papers develop this technique for higher-order derivatives and 
powers of the dependent variables illustrated with several applicable 
examples. 

The purpose of this paper will be to present the general method as 
developed by Baycura and to extend this method to varied and more com- 
plicated equations of non-linear systems. The latter will be accomplished 
by developing the transforms for the product of derivatives of the depen- 
dent variable and the dependent variable raised to the ae power, and 
derivatives raised to the A power. The validity of the transforms will 
then be shown in several applicable examples. In the examples a 
computer method designed to find the inverse Laplace Transform will be 


introduced and explained. 


II, LAPLACE TRANSFORM CONCEPT 


A. THE METHOD OF INTEGRATING BY PARTS 

The following generally summarizes Baycura's method of utilizing 
Laplace Transforms on non-linear systems using the method of integration 
by parts. Considering the product of two time-dependent functions, x(t) 
and y(t), which are to be integrated using integration by parts n - times, 


obtain the general formula: 


Cc oO ca $) (ee) 
| xydt = x } vat =e Gliteaiteet me iy VOtGtat o= coisa: 
(-1) x 2) TC | yath? ) Gite... (2-1) 
0 
fi = lc one 


where x M2) indicates the (n - 1th derivative of x, 
Applying this general formula to the Laplace Integral, equation (2-1) 
can be used to find the series form of the Laplace Transform of a function 


x(t) allowing y(t) = e~St. This integration gives the series form of: 


oe) oO oO oe) 
|seo = x | emciad: — x fos dt dt+x | J fevStatatat ~... (ae 
8 fe) o) ‘o] 


Evaluating at the upper limit it is found that the value of the integral 
approaches zero as t approaches infinity. Then evaluating at the lower 


limit the integral becomes: 


oO 


‘oes x(0 x(0)  -¥(0 
|xe st . (0) + a, 20) fT (2-3) 


or 


Siete) = x(s) = 20) + Bdge BO (2-4) 


Be APPLICATION TO DERIVAGIVES 
Equation (2-4) depicts the Laplace Transform of x(t) in an infinite 
series. expansion involving the initial conditions of the function. Through 


Similiar manipulations the transforms of the derivatives may be obtained: 


S (sabe) ] xenstg, = 4 Pea) a) 
S S S 


sX(s) - x(0) (2-5) 


Ay Betty] =) 3 eSt at = #0) +P at) + 


fe) 


= 5“ X(s) - sx(0) - x(0) (2-6) 


C. APPLICATION TO NON-LINEAR FUNCTIONS 
The simplest non-linear function to investigate is x" (t), Using the 


same technique, expansion by parts yields: 


ae (ale x jpe™ dt - x| J xe7St at at ie x| [xe SF atatat ey 


(2-7) 


Taking each of the above terms individually and again expanding by parts 


obtain: 
x|xeSt at = x(0) X(s) (2-8) 
- #ai| \ieaeieet Gude = #10) 2X5) (2-9) 


mT 


SJ xe7S* dtdtdt = #0) 26) (2-10) 


OG 


The summation of these terms then yields: 
£ ix*(t)] = | (0) + 20) a5 a Pe aera ] Ms) =e (Sy (2-11) 


Through similiar computations, transforms of higher-order functions 


may be obtained. The transform of the function cubed would be: 


i (ae |x (0) + x00) + ay lee | S42) Secs) (2-12) 
Generalizing it can be stated that: 


Deo tt) | =e eae) (2-13) 


III, LAPLACE TRANSFORMS FOR NON-LINEAR FUNCTIONS 


A. TRANSFORMS FOR THE PRODUCT OF A FUNCTION AND ITS FIRST 
DERIVATIVE 


Transforming x(t) x(t) by applying integration by parts the following 
derivation is obtained: 


Loe) o9) 
£ (x(t) x(t)] = [xxeStat = x[xe7Stat - x [fxeS! atat 


te-1) 


CO 
4 | [Jeers Giwaiwat - ... 


is) 


fet: I, = x | xe *S*at 


a oe 


Knowing the transform of the term under the integral, this term then yields: 
I, = x(0) X(s) (2) 


Similiarly allow: 


ioe) 


In= -¥ | J xe7Statat = $4(0) AG) (3-3) 
and: 
Iz = # | |Jxe7Statatat = 55110) ag) (3-4) 


The summation of these terms gives: 


iissp=e > Sn Sees 
n=l 
ses) {x(0) + 20 F (0) 4 (3-5) 


Simplifying: 


£ (x(t) x(t) ] o*x(a) |- es ee 


Ss” Xie) |X(s) ~ vil 


sX(s) [sX(s) - x(0)] (3-6) 


Continuing in this fashion the Laplace Transform for X(t)x* (t) is 


derived. 
oO oO fa) 
eee sees) | ex°e7? rai = | x2e7St gt - x | Jx?eatat 
fe) 0 ie) 
apse rr cnctaci = 50. 
2) 
co CO oO 
Let: I, = &| x®e7Stat = x x [xe7Stat - x [JxeStatdt + .... 


Further let: J, D6 |x eStqt = x(0) X(s) 


S 


iz ae x | J xe St gt dt = X(0)_ X(s) 


ja = % IJ xe~Stdtdtdt = CL EG) 
Summing: 
= = X(0) _ (0) 
IL =x oy J, = x(0) X(s) | (0) 1 eee ee | 


and after simplifying: 
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I, = sx(0) X*(s) (3-8) 
Similiarly let: 


I= - ¥|Px?e7Statdt = - x |x [fxeStatat - x [ffxe7Statatat +... | 
o 


ie) Q 


Is, = x |f xe S'atat = 210) XG) 
6 


igs = se JJ xe7S* atatat =~ AO) XS) 
Q 
Summing again: 


- *(0) X(s) j= i) AD. . 


= eae = 
a ne 


= so) X"(s) (3=3) 


By similiar arguments obtain: 


ee 


ss eee 2 
= H [ffx om *t atarar — SE) (3-10) 


Combining the above integral terms the transform is then: 


£ lsafie) x(t) | oF as OG) a a | 


n=1 


S Oise isos) — sata) (3-11) 
In similiar fashion the transform for x(t) x°(t) may be obtained and 
expressed as: 


s (x(t) x°(t)] = [x8 at = 5° x%(s) [s X(s) - x(0)] (3-12) 


o 


Equations (3-6), (3-11) and (3-12) provide the Laplace Transform for 


the general form of the function, X(t) x™(t). 
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£ (R(t) x™(t)] = sD xP(s) (£ikit)] } (3-13) 
B, TRANSFORMS FOR THE PRODUCT OF A FUNCTION AND ITS SECOND 

DERIVATIVE 

Employing again the method of integration by parts to the non-linear 


function, X(t)x(t), the transform may be expressed as: 


© joe) foe) 
Cis eadsn. = | ¥xe7Stat = ¥% [xe7Stat - ¥ | fxe St atat + ... (3-14) 


Attacking each of the integrals separately as in Section III.A we let: 


fo <) 


L = elcome = (0) X(s) (3-15) 
I, = - x {J xe"Stdtdt = #10) X65) (3-16) 
(4) a x (4) 0) X(s 
Ip = x") [lJ xeStatatdt = 210) ts) (3-17) 
Therefore: 
: ‘ : (4)(0) x(5) 
£ Bett) x) = £ ty = (0) x(ey + MOMS) , A Ae) 


Collecting terms and simplifying: 
£(K(t) x(t) ] = s X(s) [Is* X(s) - sx(0) - 2(0)] (3-19) 


The derivation of the Laplace Transform of X(t) x(t) may be similiarly 


obtained. 
SC iiscat) | = |x x®e7St gt = x | x?e7St dt Lk |{xfe7S" at dt +... (3-20) 
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Knowing the transform of x*(t) again write separately for each of the 


integrals: 
= 2 Seist _ ee 2 
I, = x [x estat = (0) sx®(s) (3-21) 
_ ose ee 2_~st owe 2 
Bo= = x | x2e -'gdtdt = (0) X2(s) (3-22) 


oe (4)(q) x2 
I, = x4) [|Tx#e7St atatdt = <0) Xs} 


fe] 


(3-23) 


Summing the above terms the transform becomes: 


o{x(t) x*(t)] = Z_ In = ¥(0) sX*(s) + 340) X*(s) wi oils) + 


S s* 4 


= stxa(s) (- 2 - 2D , 0), 2) 2) (a2) 


Pee =e xs) (Ss xis) = S10) = ea) (3-25) 


The transform for X(t) x*(t) may be identically expressed as: 


sfe(t) x(t) J = | ¥xee7Stat =x] xe St at - ¥ | [x®eStatat +... (3-26) 
£K(t) x°(t)] = s° X*(s) [s* X(s) - s x(0) - x(0)] (3-27) 


Comparison of equations (3-19), (3-25) and (3-27) provides the 


general form of the transform for these functions, 


Dice by = Ss xs) le” X(s) — sx(0) - x(0) | 


s" X"(s) {£ (k(t) ]}} (3-28) 


Ihe 


C, TRANSFORM OF THE PRODUCT OF POWERS OF A FUNCTION AND ITS 
DERIVATIVES 


Combining the results of the two foregoing sections the general form 


h power and 


of the transform of the product of a function raised to the nt 
its various derivatives can be derived. From the derived transforms as 
expressed in equations (3-13) and (3-28): 

Six(tied(t)] =sstixs(s) (2 ieee 

£(X(t) x7 (t)] = sUX"(s) (£ 1K(t) 1} 
The Laplace Transform of the product of a function raised to the nth power 


and its mth derivative is therefore: 


sfx) (ty xP(t)] = s™XM(s) Ce [xy] } (3-29) 


D, THE TRANSFORM OF THE FIRST DERIVATIVE RAISED TO THE nth POWER 
From Section II-B, equation (2-5) gives the transform of x(t). 
£ [x(t)] = s X(s) - x(0) 
Proceeding with the first derivative squared, Ex (fee and expressing 
the transform using integration by parts, the Laplace Transform becomes: 
co Pe co 


5 {(x)?] = | est = x |x e~St gt = ||[Eoue ch dt +... (3-30) 


‘e) 6 0 


Knowing the transform of x(t) the terms of the series may be expressed 


separately. 
I, = x) xe Stadt = x(0)[sX(s) - x(0)] eon 
lp = -x [fxenstatdt = 2 fs x(s) - x(0)] (3-32) 


(9) 
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I; = x | [J xe7S" dt dtdt = ae is X(s) = x(0) | (3-33) 


PES) Ste 


| 
Ms 


[sX(s) - x(0)][x(0) + + #0) ee 


[sif(s) - x(0)]s* [oa + +a + 2 nine 


Il 


s* [sX(s) ~ x(0)] [x(s) - 2) 


Sie wey One (3-34) 


Using similiar arguments the transform of Pee) * may be derived. 


oo co & 
£ x(t) P} = [xPe7Stat = x[e7Stat-x [te Statdt +... (3-35) 
) ° 8 


Knowing the transform of [x(t)]* from equation (3-34) let: 


I, = x|xeStat = x(0) s [s X(s) - x(0)]* (3-36) 
I,=- xf x?e at dt = #0) s {fs X(s) - x(0)]* (3-37) 
I, = | [J Pe7St at at ae = s iseaey = (Onn (3-38) 


Summing the above terms: 


£{R@)P}= 2 +k + b+ 


= s [sX(s) ~ x(0) J? |x(0) + + Au +2] 


LS 


=s"|s ci - x(o))? [- 8 4 + HW, 2), ve 


= 67 (eds) =a) (3-39) 


From equations (3-34) and (3-39) the transform of [x(t)]" may be obtained. 


£ {(x(t)]2}=s®7-! (sx(s) - x(0)}2 


=s"-l (etki (3-40) 
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IV. APPLICATION OF LAPLACE TRANSFORMS TO NON-LINEAR SYSTEMS 


A, EXAMPLE OF A FIRST-ORDER NON-LINEAR SYSTEM, EQUATION 1 

The first non-linear system to which the derived Laplace Transforms 
are to be applied is one which has previously been considered by Baycura 
[Ref. 3] and Brady [Ref. 4] in the initial work that has been accomplished 
in this area. The differential equation to be considered is of the form: 

x +ax + bx*® = 0 

The reasons for uSing this particular equation, aside from providing 
some measure of continuity in previous work, are thatit represents a 
stable system and is by far one of the simplest non-linear equations 
with which to work. 

In this initial analysis the author will bring to light several concepts 
and techniques that utilize Laplace Transforms to derive solutions to 
various equations and then compare these to solutions obtained from more 
“conventional” techniques. The general concept will be to express the 
transform of the unknown variable in a series of terms ins, the solution 
to which may easily be obtained on the digital computer. The two solu- 
tions used for comparison will have been obtained from an analytical 
method and by computer using Runge-Kutta fourth-order integration methods 
using Adams~-Moulton Predictor Corrector with error check, 

Utilizing, therefore, the equation: 

x + ax + bx? = 0 (4-1) 


and applying the derived Laplace Transforms we obtain: 


At 7 


bsX*(s) + (s + a)X(s) -x(0) = 0 (4-2) 
Letting x(0) =c then: 
bsX?(s). + (s.+a)X(s) = c = 0 (4-3) 


Using the quadratic formula, solve for X(s). 





2! 3 | _ 


a 
sta Abcs 
elie ala lal eal? (4-4) 
The radical may be expanded in a binomial series of the form of: 
sane (n-1) (n- 3 
(uj? 214 np ee aqui nee) + (4-5) 


Letting a=2, b=1, and x(0) =c =1 the series expression for X(s) is: 


Xe) -2- 8 ee eee eee eee 
s+2  (s+2)® “ (s+2)® (s+2)’ (s+2)® (st2) (s+2)> 9 
(4-6) 


To obtain x(t) the inverse transform of the individual terms may be obtain- 
ed and summed. An alternate and by far a more convenient and simpler 
method was, however, used in which X(s) was expressed as a ratio of 


polynomials ins, X(s) then took the form: 


ges ee + bs" hey es ee 


X(s) = 
Ci Pee Ss SS. LS Ss 
oD ctea GOT) eyo GS oat eee 


A devised computer technique was then used to obtain x(t), This method 
to find the inverse Laplace Transform by digital computer is thoroughly 
described in Appendix A, 

Using these computer methods a solution to the equation was easily 
obtained using any desired number of terms of the series. As a result, 
solutions to the equation were examined terminating with the second to 


the eighth term in the series, 
18 


several sets of solutions —s therefore been derived and are tabula:- 
ed in Table IV-1. These are an exact solution which is known and obtain- 
ed through analytical methods [Ref. 5], a computer solution using Runga- 
Kutta integration methods, and several computer solutions using Laplace 
Transforms but using varied numbers of terms in the truncated series. The 
Runga-Kutta solution is presented in order to show its accuracy Since this 


solution will be the basis for comparison in subsequent examples, 


TABLE IV-1 
SOLUTION TO EQUATION 1 


EXACT RUNGE-KUTTA APPROX, APPROX, APPROX, 
TIME (t SOLUTION SOLUTION 4 TERMS 6 TERMS 8 TERMS 
0.0 Pou 1.000 1,000 1,000 1.000 
ey OR SS DR SWAS ORS 7) eo 0.582 
0.4 0.352 OMsoZ 0.367 wore, 07370 
0.6 Oi223 Om223 0,247 0.249 O25 2 
IeeO 0.0950 0.0944 Ovo Opt 0.140 
1.4 0.0410 0.0414 0.0841 a esl) Loyd] OeUS sg 
Lucas: 0.0184 0.0183 0.0574 -0.464 0.0791 
de 2 0400823 0.00822 0.0370 Se A Sll Or033 9 
Zee 0.00368 0.00368 O02 16 -2.434 0.189 
Cra) 0.00165 0. 00265 Oya Ass =o oo O24! 
3.4 0.000734 0, OOOFE 5 0.00686 -5,094 1.402 
oye 0,000333 0.000327 0.00416 -7.024 3.246 
4.2 0.000153 0.000145 Oe00Ss0, —-8.028 6.407 
4.6 0.0000667 0.0000643 0.00334 -8,482 Ale 10)10.6) 
oO 0.0000300 0.0000283 OOO SNE™ —=8:5367 16,817 
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In comparing these various solutions, itis noted that the approximate 
transform solution conforms quite favorably with the exact solution using 
up through five terms of the series of equation (4-6). Truncating beyond 
five terms, the series has a tendency to diverge, the direction and mag- 
nitude of divergence depending on the number of terms beyond five that 
are utilized, Since this series should converge for all values of , it 


can be concluded that the original derivation of the transform: 


ty bc 2) + a) + oc +... is invalid for large values of t and 
is the direct cause of this divergence. For small values of time, t, 
however, all solutions closely approximate the true solution. 
B, EXAMPLE OF A NON-LINEAR EQUATION WITH SQUARE OF DERIVATIVE, 

EQUATION 2 

Considered now is the Froude Equation in its general form: 

K¥ +20 (l+ax)x+x=0 (4-7) 

This equation demonstrates non-linear damping, but otherwise behaves 
in a rather linear fashion. Fora stable system select. ( =0.1,@—=05 
Then equation (4-7) becomes: 

“+ 2k + .1(%)* + x = 0 (4-8) 
Taking transforms we have: 


s°X(s)e- sx(0) - #(0) + .2sX(s) — wx(0) + . loxP ie — .2x(O)eeeee 


+ .1[x(0)]°s + X(s) = 0. (4-9) 


Collecting terms and letting x(0) = 1 and x(0) =0, 


,1s°x”(s) + [8s* + .2s + 1]X(s) - (.9s + .2) =0 (4-10) 
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Now solve for X(s) using the quadratic formula and arrange ina 


Suitable form: 


%x(s) = 8s a s+] f-14[1+ As 9s + ,2) F} (4-11) 


a (.8s° +.2s5 + 1)* 
Expanding the term within the radical in a binomial series and after 


collecting terms and rearranging, X(s) can be expressed in a series. 


3 2 6 3 
wie) = eos +) ls G9sit 92) A H02S (SisehZy 


(.8s°+ .2s +1) (.8s7+.2s+1)® © (.8s7+ .2s+1)> 9 °° 


(4-12) 

The inverse Laplace Transform, x(t), is now obtained using the 
method depicted in Appendix A. The solution, compared with the equa- 
tion solution using a Runge-Kutta integration method, is displayed in 
Table [v.-”. 

The Runge-Kutta (true) solution and the approximate solution with 
values corresponding to those displayed in Table IV-2 are depicted in 
graphical form in Figure 1. The true solution is a damped sinusoid. The 
approximate solution in this case is fairly accurate for small values of 
t (up to 3 seconds), This conforms with the analysis previously given in 
section IV-A. 

Introducing appropriate initial conditions at periodic intervals of 
approximately three seconds, several piecewise but continuing solutions 
were obtained, This method provided a continuous solution but one that 
could be obtained over short periods of time. Figure 1 shows that a very 


accurate solution could be obtained in this manner. 


Zi 





TABLE IV-2 


SOLUTION TO EQUATION 2 


RUNGE-KUTTA TRANS FORM 


TIME SOLUTION SOLUTION 
Dent) 1 . 000) 1.000 
0.4 0.922 On 890 
ORs: Os 09 On735 
Th 0. 39m 0.466 
EO Oa 0kais "=. tau 
bie ili oe | BO 2) S18 
2.4 =O07o990 BUT OS 
ips sa al = Ola 
3.2 slater’: —Uroun 
3.6 -0,743 Ui. 
4.0 =('7oo6 -0 493 
4.4 cs I ee Jits' = neta 
4.8 ar Oy 0.054 
ree O22 OR AUS 
9.6 0.405 0.400 
oe, Dae Oey 20 
6.4 0,946 roo 
6.8 0.487 ORR oe is 
nee Wows 0.407 
ES OPES Z 0.249 
oO Scie 0.061 
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In contrast to this method Figure 2 depicts a true and approximate 
solution, the latter being obtained in a continuous fashion from zero to 
eight seconds. It is noted that prior to three seconds of time the solution 
is fairly accurate but subsequent to this time begins to diverge and 
becomes slightly out of phase. The divergence is noted to increase with 
tome. 

C. EXAMPLE OF THE PRODUCT OF THE VARIABLE AND ITS DERIVATIVE, 

EQUATION 3 

We next consider the van der Pol Equation which describes various 
physical situations, the more noted representing an electronic oscillator. 
The general expression for this equation is: 

¥ - a (1 -x*)x + kx = U 
Where U represents the input and @ and k are contants. Consider the 
system with zero input and with k=1, then the equation becomes: 

%®-axt+axx +x = 0 (4-12) 

Apply derived transforms, 


s* X(s) - sx(0) - x(0) -asX(s) +a x(0) +a s°X7(s) - a 8?x(0)X*(s) 
+ x{(sy = 6 (4-13) 


After collecting terms the expression becomes: 


a s°X3(s) - ws*x(0)X*s + [s® - as +1]X(s)-[sx(0)- ax(0) + X(0)]= 0 
(4-14) 


Solving this cubic equation for X(s) is no easy task; however, if we 
assume that the initial conditions for the higher derivative terms are 


negligible and equally zero then: 


Zo 


cor ees ee } (4-15) 


S S 
Substituting for X(s) in the cubed term yields: 


asx(0)X*(s) + [s* - was + 1] X(s) - (sx(0) - ax(0) + x(0)] = 0 (4-16) 


x(0) =a 
x(0) =b 
(0) - ax(0) = 


Substituting these values in the equation for X(s) a quadratic equation 
is obtained: 
absX*(s) + [s* - s+1]X(s) - fastc]=0. (4-17) 


Using the quadratic formula now solve for X(s). 


eee 2 2 : 
(s°-as + 1) +((s* -as + 1)% + (4bas) (as +c) F 
Zabs 


X(s) 


a —_ 

_s as +1 {-1+[1 + CBee) sire }} 
2abs = (s® - as + 1)* (4-18) 

Expanding the term under the radical in a binomial series and collecting 


terms yields: 


X(s) = Scene _ (bas) {as + c)* F 2(bas)*(as + c)* 
(s -as+1) (s° -as +1)? (s? - as +1) eee 


(4-19) 
Now the original equation, (4-11), represents a system having variable 
damping depending on the displacement x and hence the coefficient of the 


x term. Allowing: 


a= l 
x(0)=a=1 

x(0) =b=1 

x - wx(0) = 0 
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equation (4-19) becomes: 


S ae er 


aS ~ He? — 5 4), = aa eie me See eee 


(4-20) 


Again solve for x(t) using the devised computer methods. The results are 
Compared with a Runge-Kutta solution in Table IV-3. For asmall, as in 
our example, a number of cycles of oscillation are required before the 
steady state is achieved. The true solution is shown over a period of one 
cycle of the sinusoidal solution but during this period the initial growth 
can be detected commensurate with the slight increase in amplitude with 
appropriate increase in time. 

It should be noted at this time that the derived series representation 
for X(s) as depicted by equation (4-20) is an unstable system. This fact 
is borne out by the approximate solution as shown in Table IV-3. It can 
be generally stated then that the technique adopted here for finding 
Laplace Transforms is not applicable to systems, such as represented by 
the van der Pol Equation, which change from a condition of instability to 
one of stability. The method being investigated will not show the true 
system but will cause the solution to diverge without bound with 


increasing time, 


lua 


TABLE IV-3 
SOLUTION TO EQUATION 3 


RUNGE-KUTTA 
SOLUTION 


Le O.0 
LG 
iOS 
ldo: 
1.349 
ys ie 
0F9 072 
0.421 
= Oud. 
sls Zyl 
— ealial) 
a RS IOS 
-1,744 
= iS 
= a9 
=O 7601 
Oat) 
Ve37 
eos 


1,994 


Laes7 


TRANS FORM 


SOLUTION 


Lce,0,G 
1.174 
ive. 2- 743 
geet 
LieeOitQ 
0.071 
-2.049 
Sage 7A 
~ Gin 23 
mil SyeeeZ 3 
-16.041 
sl Sm SS 
-4.184 
MM eWanoyeat 
S02. oon 
STO, 
1S Fao, 
PADS IS: 
245.178 


PSO Hille 


148.605 


V. CONCLUSION 

In the foregoing analysis a logical method was presented by which 
the Laplace integral was expanded into an infinite series and after inte- 
grating by parts the infinite series developed could then be applied to 
non-linear systems. It was shown that this technique is applicable to 
most systems or differential equations, the exceptions being systems that 
are unstable, The limitations restricting the use of this method arise 
primarily from the complexity of the equation and the inherent difficulty 
in solving for X(s), the transform of the dependent variable x(t). The 
procedure used in this paper was to restrict the expression for X(s) to 
a quadratic equation and to further express the transform in a converging 
series uSing a binomial series expansion. 

It was observed that in obtaining the inverse of the transform a close 
approximation to the true equation/system solution could be obtained. 
This approximation was found to be valid for small values (several sec- 
onds) of solution time but for larger values of time a slow digression from 
the true solution was obtained, This anomaly was attributed to the errors 
introduced through the use of the transform expressed ina series of terms 
consisting of the initial conditions of the unknown variable and taking the 


form of: 


x(s) = #0 +a +80, 


* &@ @ 


In the attempt to find a reliable yet uncomplicated method to produce 
the inverse of the Laplace Transform and hence the solution to the system, 
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a computer program was developed whereby these painstaking computa- 


tions could then be readily handled on the digital computer, This pro- 
cedure permitted a more detailed analysis of the problem otherwise 
unobtainable through ordinary analytical methods, 

Finally, since the approximate solution was found to be valid for 
only a small interval of time a more meaningful comparison with the true 
solution was necessary. Through the introduction of appropriate initial 
conditions at periodic intervals, several piecewise yet continuing solu- 
tions to an example equation were obtained. This was easily executed 
with the use of the associated digital program, The solution obtained 


in this fashion proved very accurate and usable for all time. 


30 


APPENDIX A 
COMPUTATION OF INVERSE LAPLACE TRANSFORM 


The solution of example non-linear equations through application of 


derived Laplace Transforms is based primarily on computational techniques 


proposed by Ward and Strum [Ref. 6]. These techniques are based on 
expressing the system's equations in simultaneous first-order state equa- 
tions such as; 

X(t) = Ax(t) + But) (A-1) 
where: x(t) is the matrix of state variables 

and u(t) is the matrix of system inputs. 

Incorporating a time domain solution of these state equations, x(t), 


computed at any particular time, t, may be then expressed as: 
t 
x(t) = At x(0) + AL | EAE Bu (9 dr (A-2) 


A computer program is presented by the authors to solve for x. 


With the transform of xX given as: 


bs" + bes” a cars oa 


ni n-l n-2 
S + a, 5s + as5S + eeer + an 


is) = (A-3) 


the problem then is to transform equation (A-3) into the form of equation 
-2). 

By cross multiplying 
(s” ac tsly seek + ule +..e. t aay X(s) = b, ‘et =. ecu. ec. 


(A-4) 


ol 


The right side is then written in matrix form. 


go ,n-4 gus eT | ee (A=5) 


If X(s) resulted from the differential equation: 


(n) i=) (n-2) 


> le cen + aoS ss et + a,x = 0 with initial conditions: 
x (0) Pe 0) renee Pee (0) then the transform of the equation is: 

s'X(s) = een - _n@2,,(2) Es —— g°x lig 

+ a, [577 *x(s) - ahaa) - oo 3g (0)caieer, - -s° a2) (0): 

aly Pataes 62 

+ a) X(s) - s°x(0)] 

7 a_X(s) = 0 (A-6) 


Now collecting all I.C. terms and moving them to the right, the 


right side of the equation then becomes: 


= fle a Age ie 7 a,st79 phioeey , week aie, jwe(0) 
ae foes ze ay ire Beret ene See 
ge cg variate ts 
Us 
al ee (0) (A-7) 
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Rearranging again into a matrix product: 


= “gah ,n74 





Pree th | 0 0 ... ©O| | sae) 
| re: x) (9) 
ay | : 0 
ay Sis 
| @ 
; | 
a isi il fe § 
n-2 
(A-8) 
Now equations (A-5) and (A-8) may be equated so that: 
b, 1 0 ®  ... 0 x(0) 
be ay ] OS» cys ay (0) 
2 
; = ap a i ee 0 ) (0) 
(n-1) 
b_ a) 3n-2 - oe 1 x (0) (A-9) 
= 
For easier notation the above matrices may be written as: 
ay — ai) (A-10) 
or; x(0) =a™b (A-11) 


Applying these results to equation (A-2) note that in this case the 


input, the u matrix, will be zero, leaving but one term for x(t). 


x(t) = e&* x(0) (A-12) 
where Ais a matrix of constants as given in equation (A-1) and is of the 


form: 





JO 


=a =2 ° t ° * e s cao | 
] 


Knowing a in equation (A-11) x(0) may be easily obtained with a compu- 
ter utilizing a library program to obtain the inverse of a matrix and then 


multiplying the two resulting matrices, 
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APPENDIX B 
LIST OF LAPLACE TRANSFORMS 
£{x(t)] = X(s) 
Pie (Dae=ss i (8) 


s{x*(t)] = s* x*(s) 


o{x™(t)] = s™* x) 
£(x(t)] = sX(s) - x(0) 
fist) | = s” X{s) - s x{0) - x(0) 


en eee eg) ee 


x(0) -....-sSx (0) 
-xM-l 9) 
£ik(t)] = sX(s) - x(0) 


Lf (t)}" ] = sis X(s) - x(0)}? 
SC &(t)}" 1] = oo [sX(s) - x(0)]" = ee (£(x(t)}" 
sl oe (97) = sO) Cote ™ Gy) 


£{k(t) x(t)] = s X(s) ts X(s) - x(0)] 


g{k(t) x°(t)] = s*x°(s) [s X(s) - x(0)] 


s{x@) x @)] = s°X (s) { sik) 1} 
£e(t) x(t) ] = sX(s) [s* X(s) - sx(0) - x(0)] 


Sti) x*(t) |] = s*xX® s[s?X(s) - s x(0) - x(0)] 
s (x(t) x" (t)] = s” x" (s) (2) 1} 


ot (ty x (t)] = s” x%(s) ot ™ (13 
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